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Abstract 



The objective, of the present work, is to give a brief account on the 
transformation of a system of a higher order ordinary differential 
equations to another system of first order linear differential equations 
with certain applications. We are interested in the case of constant 
coefficients. We started by presenting some guide informations about 
the classification of the first order linear systems of equations. Then we 
explained how to transform a single second order linear differential 
equation to a system of two linear differential equations of the first 
order and we gave some applications on the different cases. 

We considered the transformation of a system of two second order 
linear differential equations to a system of four linear differential 
equations of the first order. This is done through the physical application 
of the coupled pendula. We determined the matrix of coefficients and 
calculated its eigen values and the corresponding eigen vectors. The 
solution of the original problem is obtained and drawn. 
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Introduction 



We will now turn our attention to solving systems of simultaneous 
homogeneous first order linear differential equations, we have to 
limit ourselves to the simplest instances: those systems of two 
equations and two unknowns only. But first, we shall have a brief 
overview and learn some notations and terminology. 
A system of n linear first order differential equations in n unknowns 
(an n xn system of linear equations) has the general form: 
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X 
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/ _ 



/ _ 



X 



CI nX i + ^12^2 ~*~ 
a 3l X l Cl 32 X2 + 
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Where the coefficients ay's, and^/'s are constants. If every 
Term Q\ r is constant zero, then the system is said to be homogeneous. 
Otherwise, it is a nonhomogeneous system if even one of the g's is 
nonzero. 

The system (*) is most often given in a shorthand format as a matrix- 
vector equation, in the form: 
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Where the matrix of coefficients, A, is called the coefficient matrix of 
the system. The vectors x' , x, and g are 



x — 



rxi 



X-. 



X — 



Lx r 



x 3 



Lx n J 



g = 



92 

9i 



-9n- 



For a homogeneous system, g is the zero vector. Hence it has the 
form: 

x — Ax 
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Transformation of a second order 
linear differential equation to a 
system of two linear differential 
equations of the first order 

We give some simple examples to explain the process 
Examples: 

(i) The mechanical vibration equation m u" + yu' + ku = F[t) is 
equivalent to the following system 

X 1 — x 2 

, _ -k y F(t) 

Xn — X\ X? "I 

-m -m w 

Note that the system would be homogeneous (respectively, 
nonhomogeneous) if the original equation is homogeneous 
(respectively, nonhomogeneous). 

(ii) The equation y'" - 2y" + 3y' - 4y = is equivalent to following 
system 

xi - X2 
X2 - X3 

X3' = 4 XI - 3 X2 + 2 X3 

This process can be easily generalized. Given an n-th order linear 
equation of the form 

CLny{n) + CLn-iy{n-l) + CLn-2y{n-2) + ... + diy" + fllj/ + QOJ/ = g(t). 

Make the substitutions: xi=y, xi -y ' , X3 = y" , ... , xn = y{n-i), and xn =y{n). 
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The first n - 1 equations follow thusly. Lastly, substitute the x's into 
the original equation to rewrite it into the n-th equation and obtain 
the system of the form: 



xi - X2 
XI - X3 
X3' - XA 



Xn-1 - Xn 



■n — -A- 1 Ji 2 JL 3 — ••• Ji n H 

^-n, &<ri 



Exercises: 

1-3 Convert each linear equation into a system of first order 
equations. 

1. y" - 4/ + 5y = 

Answers: 

1. Xl' = X2 

xi = -5xi + 4X2 

2. j/"-5y" + 9y=tcos2t 

Answers 

2.Xl' = X2 

X2' = X3 

X3' = -9x1 + 5x3 + t cos 2t 

3. y( 4 ) + 3/" - tt/' + 2ny' - 6y = 11 

Answers 

XI ' = X23. 
X2' = X3 
X3' = X4 

X4' = 6X1 - 27TX2 + 7TX3 - 3X4 +11 
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Solution of a Homogeneous System of Two Linear 
First Order Equations in Two Unknowns 



Given: 

X' — AY* 



We first find the two eigenvalues, r, and their respective 
corresponding 

eigenvectors, k, of the coefficient matrix A. Depending on the eigen 
values and eigenvectors, the general solution is written. Three 
different cases arise 

Case I. Two distinct real eigenvalues ri and r2 of the matrix A: 
The solution is 

X = CiM 1 "^ + c 2 k 2 e r ^ 

Example '■ 




This coefficient matrix has the eigenvalues r= -1 and 6. And they 
each respectively has an eigenvector 




Therefore, a general solution of this system of differential equations 
is: 

X = c i[_ 1 1 ] e " t + c 2^]e 6t 

Case II. Two complex conjugate eigenvalues A ± where A + \ii has as 
an eigen vector k = a + bi: 
The solution is 

X = C!e At (<^ cos(nt) - b sin(nt) ) + c 2 e Xt (a sin(|it) + b cos(^it) ) 
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Example '• 



\2 -51 
X = ll -2\ X - 



The characteristic equation is r 2 + 1 = 0, giving eigenvalues r = ±i. 
That is, A = and \i =1. 

Take the first (the one with positive imaginary part) eigenvalue r = z, 
and find one of its eigenvectors: 

Solving the first equation of the system: (2 - /) xi - 5x2 = 0, we get 
the relation (2 - /) xi = 5x2. Hence, 



Therefore, the general solution is 



X = 



cie 0t ([5] cos(t) - [0J sin(t) ) + c 2 e 0t ( 
+ [Oj cos(t) ) 



sin(t) 



5 cos(t) 
2 cos(t) + sin(t) 



+ c 2 ( 



5sin(t) 
2sin(t) - cos(t) 



Case III. A repeated real eigenvalue r: 
There are two cases 

(Y) When two linearly independent eigenvectors 
exist - 



= Cl k ie rt + c 2 k 2 e rt 



11 



Example '• 




The eigenvalue is r = 2 (repeated). There are 2 sets of linearly 
independent eigenvectors, which could be represented by any 2 
nonzero vectors that are not constant multiples of each other. For 
example: 

* = 0. = 0- 

Therefore, a general solution is 

X = d [J] e 2t + c 2 g] e 2t 
(ii) When only one linearly independent eigenvector exist - 

X = Cl k ie rt + c 2 (kte rt + T]e rt ) 
Note '■ Solve the system [A - r l)t] = kto find the vector t). 



Example '■ 



; = [4 -7]* < *°> = [ i \ 



The eigenvalue is r = -3 (repeated). The corresponding system is 



Both equations of the system are 4xi- 4x2 = 0, we get the same 
relationxi = X2. Hence, there is only one linearly independent 
eigenvector : 

* = a- 

Next, We solve for 17: 
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[I ill HI] 



It has solution in the form 
rl + T] 2 



Choose T]2 = 0, we get 
1 



4 

OJ 



A general solution is, therefore, 




Apply the initial values to find Ci = 1 and Ci = -12. The 
particular solution is 



X 



-0 



-3t 



- 12 




-3t 



+ 







-3t - 



-12t-2 
-12t+U 



-3t 
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Transformation of a system of second 
order linear differential equations to 
a system of four linear differential 
equations of the first order 

We explain this process through the physical application of the 
coupled pendula given here after 



14 



The coupled pendula 



Two identical simple pendula oscillate in the plane as shown in 
Figure 1. Both pendula consist of light rods of length ■£ = 10 and are 
suspended from the same ceiling a distance L = 15 apart, with equal 
masses m = 1 attached to their ends. The angles the pendula make to 
the downward vertical are x and 2 , and theyare coupled through 
the spring shown which has stiffness coefficient k = l.The spring has 
unstretched length L = 15. 

We assume that the acceleration due to gravity g = 10. 

(a) Assuming that the oscillations of the spring remain small in 
amplitude, so that 1 0x1 «land|0 2 | «1. 



Application of Newton's second law and Hooke's law, shows 
that the coupled pendula system gives rise to the system of 
differential equations 



It is usual in the literature to put this system of second order 
differential equations in the matrix - vector form 




-20! + 9 2 
X - 26 2 



(1) 
(2) 




= Ay 



(3) 



dt 2 
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Where 



-J 



and y = 




Here Sis known as the stiffness matrix and y is the vector of 
unknown angles for each of the pendula shown in Figure 7.1. 

(b) By looking for a solution of the form 



for a constant vector c, show that solving the system of differential 
equations (3) reduces to solving the eigenvalue problem 



(c) Solve the eigenvalue problem (4) in part (b) above, stating clearly 
the eigenvalues and associated eigenvectors. 

(d) Hence enumerate the possible modes of oscillation of the masses 
corresponding to each eigenvalue-eigenvector pair. 

(e) Finally, write down the general solution of the system of 
equations (3) 



y(t) = ce 



At 



(B - A 2 I)c = 



(4) 



16 




\ 




\ 



m=l 



FIGURE 1. Simple coupled, peridula system. 



Transformation to a linear system of the first 
order 

In the following, we study this problem by transforming the above 
system of two second order differential equations to an equivalent 
system of four linear differential equations of the first order as 
follows 

Solution: 

We reconsider the system of second order equations : 

d\ — — 2.0i + O2 
2 = Q 1 — 26 2 

And we define 
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This leads to => 

— — "ZOi + 2 



64 — 6 1 — 26 2 

Thus, we have produced a system of four equations and four 
unknowns, 6 1 ,6 2 , 3 ,6 4 



We consider that we have a vector^ 



Where 



x = 



0i 
4 



X = AX 



- 


1 


0- 








1 


-2 


1 





- 1 


-2 


0- 



To solve this system we find eigen values and eigenvectors of the 
matrix A. The four eigen values are: 

^1,2 — ±i , ^-3,4 — ±V3i 
Mfe usee/ the MATLAB software to find the eigenvectors as 



W-l = U-l + iV-L 



-0.5 
-0.5 
-0.5i 
-0.5iJ 



and w 2 — u 2 + v 2 — 



r-.3536t 

0.3536i 
0.6124 
L- 0.6124J 



p = [v 1 u x v 2 u 2 ] 
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r -0.5 -0.3536 

-0.5 0.3536 

-0.5 0.6124 

1—0.5 -0.6124J 



Then 

r -1 -1 ] 

-1 = -1 -1 

p ~ -1.4140 1.4140 

0.8165 -0.8165- 



P~ X AP = 






-1 








1 




















-V3 








V3 






The solution subject to the initial condition 

x(0) = x Q 



is given by 



x(t) = P 



'cost —sint 

sint cost 

cosV3t —sinV3t 

. smV3t cosV3t . 
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From which, we get the final result as 



x(t) = 

-O.Scost + 0.4cosV3t -O.Scost - 0.4cosV3t 

O.Scost - 0.4cosV3t O.Scost + 0.4cosV3t 

O.Ssint — 0.6sinV3t O.Ssint + 0.6sinV3t 

-— O.Ssint + 0.6s£nV3t — O.Ssint — 0.6s£nV3t 



-O.Ssint + 0.4smV3t -O.Ssint - OAsiny[?>t 
O.Ssint + 0.4sinV3t O.Ssint - 0.4sinV3t 

— O.Scost + 0.6cosV3t — O.Scost — 0.6cosV3t 

O.Scost + 0.6cosV3t — O.Scost — 0.6cosV3t - 



The vector x represents the initial values of the functions 6-1,62, 6 3 , 6 4 
and the product of the matrices gives the values of these functions at 
any instant of time t. 

We notice that the third and the forth rows are respectively the 
derivatives of the first and the second rows and this result indicates 
the exactness of our solutions since from the beginning we assumed 
that 6> 3 = 6\ and 6 4 = 6' 2 . 

The following four functions , appearing as a first row of the square 
matrix, determine the solution for any initial values of the functions 
9 lt 2 , 3 , 4 between t=0.0 and t=20 

Fi(t)= -O.Scost + 0AcosV3t, 
F2(t)= -O.Scost - 0Acos43t, 
F3(t)= -O.Ssint + 0AsinV3t 
and 

F4(t)= -O.Ssint - 0AsinV3t. 
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ll 1 1 1 

0.5- 
Fl(t) 0- 
-0.5- 

_ 1 \ 1 1 1 

5 10 15 2C 

t 

Figure 2: The function Fl(t) 
ii 1 1 1 

0.5" 
F2(t) 0- 
-05- 

_ l \ 1 1 1 

5 10 15 2 

t 

Figure 3 : The function F2(t) 
ii 1 1 1 

05- 
F3(t) 0- 
-05- 

_ 1 \ 1 1 1 

5 10 15 2 

t 



Figure 4: The function F3(t) 




Figure 5 : The function F4(t) 



21 



Conclusions 



The main results of the present work are 

(1) presenting how to transform a higher order linear 
ordinary differential equation to a system of first order 
linear differential equations. 

(2) giving certain applications to make clear the step (1). 

(3) Presenting how to replace a system of second order 
ordinary differential equations by a suitable system of the 
first order. 

(4) Giving a physical application of the coupled pendula as 
an illustration of the method and presenting two methods 

for the solution. 

(5) giving the graphs of the functions determining the 
solution of the coupled pendula for any given initial 
situations. 
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